t = 0:1:100;
var_r2 = zeros(100,1);
var_r2(1) = 0.1;
for i = 2:101
    var_r2(i) = varR2(zeros(20,1),zeros(20,1),0.2,[t(i-1),t(i)],0.01);
end
plot(t,var_r2);
xlabel('time(s)');
ylabel('var(R^2)');
